home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b122.asm < prev    next >
Assembly Source File  |  1992-04-28  |  4KB  |  76 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.22     NxN  NxN  Matrix Multiply  
  8. ;The matrix multiplications are for square NxN matrices. All the elements  are stored in "row major" for-
  9. ;mat. i.e. for the array A:  
  10. ;   a(1,1) ... a(1,N) 
  11. ;                                          a(N,1) ... a(N,N) 
  12. ;
  13. ;the elements are stored:  
  14. ;a(1,1), a(1,2), ..., a(1,N), a(2,1), a(2,2), ..., a(2,N), ... 
  15. ;
  16. ;The following code implements C=AB where A and B are square matrices.  
  17. ;
  18. ;                                      DSP56000 IMPLEMENTATION 
  19. ;                                                           Program     ICycles
  20. ;                                                            Words 
  21. ; 
  22. ; move   #mat_a,r0                 ;point to A             1        1 
  23. ; move   #mat_b,r4                 ;point to B             1        1 
  24. ; move   #mat_c,r6                 ;output mat C           1        1 
  25. ; move   #N,n0                     ;array size             1        1 
  26. ; move   n0,n5                                             1        1 
  27. ; 
  28. ; do     #N,_rows                  ;do rows                2        3 
  29. ; do     #N,_cols                  ;do columns             2        3 
  30. ; move   r0,r1                     ;copy start of row A    1        1 
  31. ; move   r4,r5                     ;copy start of col B    1        1 
  32. ; clr    a                         ;clear sum and pipe     1        1 
  33. ; move             x:(r1)+,x0 y:(r5)+n5,y0                 1        1 
  34. ; rep    #N-1                      ;sum                    1        1 
  35. ; mac    x0,y0,a   x:(r1)+,x0 y:(r5)+n5,y0                 1        2 
  36. ; macr   x0,y0,a   (r4)+           ;finish, next column B  1        1 
  37. ; move   a,y:(r6)+                 ;save output            1        1 
  38. ;_ecols 
  39. ;; move   (r0)+n0                   ;next row A             1        1 
  40. ; move   #mat_b,r4                 ;first element B        1        1 
  41. ;_erows 
  42. ;                                                        -----    ----- 
  43. ;                                                         19          
  44. ;                               ((8+(N-1))N+5)N+8 = N3 +7*N2 +5N+8  ? 
  45. ;At a DSP56000/1 clock speed of 20.5 MHz, a [10x10][10x10] can be computed in .1715 ms.
  46. ;
  47. ;                                             DSP96002 IMPLEMENTATION 
  48. ; 
  49. ;                                                            Program    ICycles 
  50. ;                                                            Words 
  51. ; 
  52.    move   #mat_a,r0                 ;point to A              1        1 
  53.    move   #mat_c,r6                 ;output mat C            1        1 
  54.    move   #N,n0                     ;array size              1        1 
  55.    move   n0,n5                                              1        1 
  56.  
  57.    do     #N,_rows                                      ;    2        3 
  58.    move   #mat_b,r4                 ;point to B              1        1 
  59.    move   r0,r1                     ;copy start of row       1        1 
  60.    do     #N,_cols                                      ;    2        3 
  61.    move                         r4,r5                   ;    1        1 
  62.    fclr   d0                    (r4)+                   ;    1        1 
  63.    fclr   d1                    x:(r1)+,d4.s y:(r5)+n5,d5.s ;1        1 
  64.    rep    #N                                                ;1        2 
  65.    fmpy   d4,d5,d1 fadd.s d1,d0 x:(r1)+,d4.s y:(r5)+n5,d5.s ;1        1 
  66.                    fadd.s d1,d0 r0,r1                       ;1        1 
  67.    move                                     d0.s,y:(r6)+    ;1        1 
  68. _cols 
  69.    move                         (r0)+n0                     ;1        1 
  70. _rows 
  71. ;                                                         ; -----    ----- 
  72. ;                                             Totals:       19          
  73.                                                                        
  74. ;                                                ((N+7)N+6)N+7 = N3 +7*N2 +6N+7 ?
  75. ; At a DSP96002 clock speed of 26.66 MHz, a [10x10][10x10] can be computed in .1325 ms.
  76.